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(Dated: October 3, 2000) 

Single-dot spectroscopy is now able to resolve the energies of excitons, multi-excitons, and 
charging of semiconductor quantum dots with < 1 meV resolution. We discuss the physical 
content of these energies and show how they can be calculated via Quantum Monte Carlo (QMC) 
and Configuration Interaction (CI) methods. The spectroscopic energies have three pieces: (i) 
a "perturbative part" reflecting carrier-carrier direct and exchange Coulomb energies obtained 
from fixed single-particle orbitals, (ii) a "self-consistency correction" when the single particle 
orbitals are allowed to adjust to the presence of carrier-carrier interaction, and (iii) a "correlation 
correction." We first apply the QMC and CI methods to a modle single-particle Hamiltonian: 
a spherical dot with a finite barrier and single-band effective mass. This allows us to test the 
convergence of the CI and to establish the relative importance of the three terms (i) - (iii) above. 
Next, we apply the CI method to a realistic single-particle Hamiltonian for a CdSe dot, including 
via a pseudopotential description the atomistic features, multi-band coupling, spin-orbit effects, 
and surface passivation. We include all bound states (up to 40,000 Slater determinants) in the 
CI expansion. Our study shows that: (1) typical exciton transition energies, which are ~ 1 eV, 
can be calculated to better than 95% by perturbation theory, with only a ~ 2 meV correlation 
correction; (2) typical electron addition energies are ~ 40 meV, of which correlation contributes 
very little (~ 1 meV); (3) typical biexciton binding energies are positive and ~ 10 meV and 
almost entirely due to correlation energy, and exciton addition energies are ~ 30 meV with 
nearly all contribution due to correlation; (4) while QMC is currently limited to a single-band 
effective mass Hamiltonian, CI may be used with much more realistic models, which capture 
the correct symmetries and electronic structure of the dots, leading to qualitatively different 
predictions than effective mass models; and (5) and CI gives excited state energies necessary to 
identify some of the peaks that appear in single-dot photoluminescence spectra. 



I. INTRODUCTION: THE PHYSICAL 
CONTENT OF EXCITON, MULTIEXCITON, 
AND CHARGING ENERGIES IN DOTS 

Small semiconductor dots, such as semiconductor em- 
bedded Stranski-Krastanow (SK) dots or "free-standing" 
colloidal dots, are engineered anaLi-studied for their 
optical and transport properties SMa Measurements 



Only differences in energy are accessible to experiment. 
We focus on four physical quantities: 

(i) Exciton energies. The exciton transition energy 

(ii) 

E x is the difference in total energy of a dot having 
as a dominant configuration an electron in level and a 
hole in level hj and a dot in the ground state. 
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ajjound,— quantities such 

„, JMM and charging 

l 15 H 16 l Advanced experimental techniques, 
such as single-dot spectroscopy, are able to resolve to 
such energies to < 1 meV resolution. This article dis- 
cusses the physical content of such measured quantities 
in terms of the mean-field (direct and exchange) Coulomb 
energies, which are relatively simple to model, and corre- 
lation energies, which we calculate by two leading meth- 
ods in the field — Quantum Monte Carlo (QMC) and 
Configuration-Interaction (CI). 

Let us consider a quantum dot with M holes in the 
valence band and N electrons in the conduction band. 
The total energy of the dot is Em,n(o:), where a is a 
quantum number that identifies the state of the system. 



Typical excitonic transition emesgies in III-V or II- VI 
dots, measured experimentallyQ'tH by photoluminescence 
(PL) or by absorption, are 1-3 eV. The exciton binding 
energy Ax is the difference between the total energy of 
a system consisting of two infinitely separated identical 
dots, one with a hole in ho and the other with an elec- 
tron in eo, and the total energy of a quantum dot with 
an exciton: 



Ax — £q,o + £0,1 — -E-1,1 — Ei 



0,0, 



(2) 



where stands for E\fi{h^^), E ,i — £o,i(^o e o)> and 
= £d,i(ft<o e o)- Typical exciton blading energies in 
III-V and II- VI dots are 10-200 meV.00 
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(it) Biexciton energies. The biexciton binding energy 
Axx is the difference between twice the exciton energy 
(or the energy of a system of two infinitely separated 
dots, each with an electron- hole pair), and the biexciton 
energy: 

^■xx — 2£i,i - £2,2 _ E 0tQ . (3) 

The biexciton binding energy is positive ("bound biex- 
citon") when the total energy of two excitons in the 
same dot is lower than the energy of the two excitons 
in two separate dots. A bound biexciton appears as a 
red-shifting of the exciton luminescence energy when a 
second exciton is present. This jjias^ seen in single-dot 
spectroscopy e. g. for Tn A s/ Gn A c Hnl 1 W^RWp.t^n bjud- 
ing energies in III-V dots are 1 - 6 mevMBBOO'El 

(Hi) Multi-exciton energies. The N-th exciton charg- 
ing energy Wn is the minimum energy needed to add to 
a dot having N — 1 electron-hole pairs (excitons) in their 
ground state one additional exciton, 

Wn — En,n — En-\,n-\- (4) 

Physically, Wn is the highest possible energy for a photon 
emitted in the transition from the lowest energy state of 
N excitons to a state with N — 1 excitons. The difference 
in successive multi-exciton charging energies is the N-th 
exciton addition energy ^^n+v 

AX) 

(5) 

— En+\,n+x + Eff—i,N—i ~ %En,n- 

(iv) Electron loading energies. The electron charging 

energy jj, N is the chemical potential needed to add an 
electron to a dot already having N — 1 electrons: 

/% = Eo,n — Eo r N-l, (6) 

whereas the electron addition energy is the difference be- 
tween two successive chemical potentials, 
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_ (e) 
-1 ~~ f-N 



Mjv-i 



= Eo^n+i + Eq ; n-i — 2E 0i n- 

Electron addition energies in .. rp lloidal dotsQ are ~ 
200 meV, whereas in SK dotsBEa they are ~ 20 meV. 

The definitions given here in Eqs. (Q)-(0) are opera- 
tional, model-independent. A central question in the field 
is how to approximate these quantities through models. 
This requires knowing how much of the energy involved in 
the precesses described by Eqs. (|l|)-(0) are due to "mean- 
field" effects, which can be modeled relatively simply, 
and how much is due to inter-particle correlation, which 
is more intricate to model. Figure |l| illustrates the steps 
required to model the electronic and optical properties 
of a quantum dot: (a) choosing a structure (including 
size, shape, composition, and strain), (b) solving a sin- 
gle particle model, and (c) treating interactions among 



Modeling the Electronic and Optical Properties of a Quantum Dot 
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FIG. 1: The three steps to modeling a quantum dot. (a) The 
structure is modeled by choosing a size, shape, and compo- 
sition profile, and determining the strain, (b) Single parti- 
cle properties of the electrons and holes are found by solving 
a Schrodinger equation for a chosen level of renormalization 
(EMA, k ■ p, tight-binding, or pseudopotential). (c) Inter- 
actions between excitations (electron-electron, electron-hole, 
hole-hole) are added to the single particle model, using ei- 
ther perturbation theory Ept; self-consistent mean-field the- 
ory, which adds the self-consistent contribution A_Esc, or full 
treatment (CI or QMC), which adds the correlation energy 

Ecorr • 



the electrons and holes. In this paper we are concerned 
with general trends in correlation in dots, so we focus 
mainly on the choice of single particle model, Fig. |l|(b), 
and treatment of interactions, Fig. |l|(c). 

As illustrated in Fig. |l](b) , the calculations of the quan- 
tities of Eqs. (§)-(0) require one to assume an under- 
lying single-particle model, which determines the single 
particle states (conduction electrons and holes). The 
single-particle model is cast as a Schrodinger equation 
with an effective single-particle potential. This potential 
contain all structural information about the system: the 
size, shape, composition, surfaces, interfaces of the dot 
system. Various levels of renormalization exist for the 
quantum dot single-particle model. The simplest is an 
effective mass ("particle-in-a-box") model, in which the 
electron and hole excitations come from single parabolic 
bands. Better approximations are the multi-band k ■ p, 
tight-binding, and pseudopotentials. 

The single-particle models do not usually contain the 
Coulomb interactions between the single-particle excita- 
tions (i. e. electron-electron, electron-hole, and hole- 
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TABLE I: Relationship between the choice of single particle 
models for quantum dots and the availability of QMC and CI 
methods to calculate correlation energy [see Fig. 0, (b) and 
(c)]. This information motivates our approach to studying 
correlation: first we test the convergence of CI against QMC 
calculations using a simple single-band EMA model, then we 
present CI calculations on a realistic multi-band pseudopo- 
tential quantum dot model to illustrate features missed by 
the simple model. 



Level of Rcnormalization Model 



CI" QMC 



All electron 



Valence only 



Active electron only 



Exact Hamiltonian no no 

Multi-Band Pseudopotcntial yes no b 
Tight-binding yes no b 



Multi-Band k ■ p 
Single-Band EMA 



yes 
yes 



no 

yes 



"While (Jl may be applied to any model, it is ottcn under con- 
verged. 

'Possible for very small clusters of less than 100 atoms (R < 10 A). 



hole excitations). Instead, these interactions must be 
added to the model, as shown in Fig. |l|(c). We classify 
the treatment of interaction among the single-particle 
excitations in-..trirpi? JeM&Ls: (i) first order Perturbation 
Theory (PT)j^&BS3nM which includes direct and ex- 
change Coulomb interactions, J and K, evaluated from 
fixed single-rparticle orbitals; (ii) Self- Consistent Mean- 
Field (MF)p!l in which the direct and exchange Coulomb 
terms are solved self-consistently [the difference between 
(ii) and (i) is called "self-consistency coroaddQn?L [ AES c ]; 

fhffd^sMfh as CBEMeBB or 
which include all many- 
body effects of interactions. The difference between the 
exact energy (iii) and the mean-field energy (ii) is called 
the "correlation correction," Ai5 corr . Thus, the energy 
for a dot with M holes and N electrons can be separated 
into three terms, 



and (iii) jM^^eJx^iLj^ihffd^—iUfh as 

QMc^k3xmkimmm)B whk 



ptot _ pPT 



a Tpcorr 
^M.N 



(8) 



which arc perturbation theory E^n > self-consistent cor- 
rections AE^ N , and the correlation correction AE^" N . 

Due to computational limitations, the methods avail- 
able to calculate correlation are dependent on which 
single-particle model is chosen (level of renormalization) . 
The computational cost for accurately calculating cor- 
relation energies increases rapidly with the number of 
electrons one needs to consider. The number of electrons 
depends on both the dot's size and on the type of renor- 
malization one uses for the Hamiltonian. As summarized 
in Table B, three levels of renormalization are pertinent: 

(a) The all electron approach, where the number of 
electrons per atom equals its atomic number. Thus, Si 
has 14 electrons per atom, and a 40 A diameter spherical 
Si dot has 1, 600 x 14 = 22, 400 electrons. This is outside 
the reach of QMC, CI, and density functional methods. 

(b) The valence-only pseudopotential approach, where 
the "core" electrons are removed as dynamic variables 



and replaced by an (often non-local) ionic potential. 
Thus, Si has 4 electrons per atom, and a 40 A diameter 
spherical Si dot has 1, 600 x 4 = 6, 400 electrons. This is 
outside the reach of density functional methods, and too 
large for QMC calculations, which axe currently limited 
to about 25 Si atoms (100 electrons).c3 Note that the all- 
valence pseudopotential approach can be further simpli- 
fied, with no additional approximations by-seaijciung for 
eigensolutions in a fixed "energy window, "ej : EJ'l3 e. g., 
near the band edges. Thus, a 40 A diameter Si dot would 
require calculating ~ 10 eigensolutions. This jteiekjnakpi 
pseudopotential calculations of dots feasible ,E3E3cac3'c2l 
and CI may be used to comput e . co i rr fl la tion energies 
from the single particle solutions .c3E3E§E3 It would be 
interesting if such folding techniques could be applied to 
QMC. 

(c) The "active-electron-only" Effective Mass Approx- 
imation (EMA) approach, where all of the "indigenous" 
core and valence electrons are eliminated (replaced by 
dielectric screening) and only additional, band-edge elec- 
trons and holes are considered. Thus a 40 A diameter 
Si dot has zero electrons. One can study added elec- 
trons and holes. This renormalization represents a severe 
approximation with respect to levels (a) and (b) above. 
Both QMC and CI methods may be readily applied to 
EMA Hamiltonians. Some improvement can be made by 
using several bands to describe the-padditional electrons 
and holes using the k -p formalism Jljo but current QMC 
methods do not treat k ■ p Hamiltonians. 

Most correlated calculations on quantum dots have 
used such a single-band effective mass model [level 
(c), above], where multi-band and inter- valley couplings 
are ignored. This particle-in-a-hqx^desjcxiDtion of the 
mean-field problem was recentlyBBMia contrasted 
with the pseudopotential solution of the problem [(b) 
above] both for "free-standing" (colloidal) dots and for 
semiconductor-embedded SK 4^ts. It.-was found that 
for "free-standing" dots (InPjla CdSeB) the effective 
mass approach can lead to energy shifts of the orderc£lcZl 
~ 500 meV; lead to reverse order of (s,p) levels;E3 miss 
more than half of the single-particle .eigenvalues in a 
0.5 eV energy range near the band edgejlliinderestimate 
the Coulomb integrals Jy [Eq. @] byB ~ 20%; and 
miss all the long-range part .of the exchange integralscij 
Kij. For pyramidal SK dotsC3 the errors are somewhat 
smaller: shifts in the energy levels for electrons and holes 
are ~ 35 meV and ~ 110 meV, respectively; energy spac- 
ings from EMA are about a factor of two too large; and 
the polarization ratio for dipole transitions along the two 
directions is 1 instead of 1.3. Such limitation in the EMA 
create a dilemma when modeling correlation as summa- 
rized in Table |. On one hand CI expansions maybe ap- 
plied to realistic single-particle models (e. g. pseudopo- 
tentials), but converge slowly with the number of con- 
figurations. On the other hand, QMC methods can give 
numerically exact answers including all correlation, but 
currently are limited to simple single-band effective mass 
models. This situation prompts us to use the following 
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strategy to study correlation effects: First, we consider 
a simplified "particle-in-a-box" single-band EMA model 
which can be treated both via QMC and CI. Our best 
CI calculations for the EMA model include all bound 
states, but neglect continuum states. Second, we consider 
a CdSe dot whose single-particle properties are described 
realistically by pseudopotentials, and the correlation is 
treaded via CI only. 

Our single-band EMA dot has been chosen to be rep- 
resentative of SK and colloidal dots. We summarize the 
properties of our model dot in Table We find that for 
the single-band model dot: 

(i) Typical exciton transition energies for our model 
dots are ~ 1 eV, and typical exciton binding energies 
are ~ 50 meV. Of this, MF gives > 95% of the binding 
energy. Correlation is only ~ 2 meV, of which QMC 
provides an accurate solution Although CI misses half 
the correlation energy, i.e. ~ 1 meV, it still captures 
~ 98% of the total binding energy. 

(ii) Typical biexciton transition energies for our dots 
are ~ 2 eV and typical biexciton binding energies are 
~ 6 meV. The biexciton binding energy from mean-field 
theory is slightly negative (unbound biexcitons), so the 
positive biexciton binding is in fact due to ~ 6 meV of 
correlation energy. QMC captures all the correlation en- 
ergy, whereas our CI captures only half (about 4 meV), 
so that the CI estimate of biexciton binding is only about 

Ealn — Eofl + (^ e c — ^ e v ) + ^ (Jv,v> 

c v v<v' 

where £j are the single-particle energies, Jij are the di- 
rect Coulomb energies, and Kij are the exchange en- 
ergies. The single-particle energies £j are often ob- 
tained from the solution of an effective single-particle 
Schrodinger equation, 

{~V 2 +V tS }il> i = e i il> i (io) 

where V c s is an effective potential. The Coulomb and 
exchange energies are given in terms of the single-particle 
wave functions ipi by: 

J e(r,r') |r-r'| 

/■ ^(r)^(r)^(rQ^(rQ , 

K i:j = / — drdr', (11) 

J e(r,r') |r-r'| 

where e is the dielectric constant of the quantum dot. 

The self-consistent contribution Efp N , given by the 
first two terms on the right hand side of Eq. (§1), arises 
from the self-consistent rearrangement of the single- 
particle wavefunction in respond to the electrostatic field, 
Eq. (pi]), generated by the excitation of electrons and 
holes. 



65% of the true value. 

(iii) Typical electron charing energies for our dots are 

(Jf* w 150 meV, relative to the dot material CBM, while 
addition energies are ~ 40 meV. Of this, correlation 
energy is very small (~ 1 meV), so mean-field or even 
perturbation theory describes dot charging and addition 
energies very well. 

For our realistic CdSe dot we find that CI can be effec- 
tively combined with accurate pseudopotential descrip- 
tion of the MF problem, thus incorporating surface ef- 
fects, hybridization, multi-band coupling. Furthermore, 
CI can calculate excited states easily, thus obtaining the 
many transitions seen experimentally, rather than only 
ground-state-to-ground-state decay calculated by con- 
ventional QMC (note, however, that-astensions of QMC 
to several excited states are possibleOH) . 

II. METHODS OF CALCULATION 

A Uncorrelated methods: perturbation theory and 
mean held methods 

The first-order perturbation energy E^ N [Eq. (|[)] can 
be written analytically as: 



Kv,v>) + ^(Jc,c> - K c ,c) - £(J„, C - K VlC ), (9) 

c<c' v,c 

B The correlated, many-particle methods 

1 Quantum Monte Carlo 

The original QMC methodi! was based on the varia- 
tional technique, a simple, yet powerful theoretical tool. 
In a variational calculation, one proposes a parameterized 
trial wavefunction ^S^(R), where A represents a set of 
variational parameters and R represents the coordinates 
of all the particles. The energy expectation value 

eW = JdR* { T X} *(R)H* { T X} (R) 
T fdRV { T X} *(R)y { T X} (R) ' 

may be minimized with respect to the variational pa- 
rameters A to give an estimate for the ground state 
energy and ground state wavefunction. This integral 
may be evaluated analytically, or Monte Carlo integra- 
tion may be used. In this simplest formulation, QMC 
is formally equivalent to the variational techniques com- 
monly applied to excitons in nanostructures.O Because 
the integral is over all electron and hole coordinates R, 
variational QMC calculations resemble classical simula- 
tions: a configuration of particle positions R undergoes 
a random walk through configuration space, using the 
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TABLE II: Measurable quantities for our single-band spherical model dot, with effective masses m e = 0.1 and mj, = 0.5, 
dielectric constant e = 12, dot-material band gap E ga _ p = 1 eV, and band offsets AE V = 200 eV and AE C = 400 meV. For each 
quantity we give the magnitude, the mean-field value, the correlation correction, and the percent of the energy recovered by 
CI expansion of using all bound states. All energies are given in meV, and electron charging and total energies are measured 
relative to the dot material CBM. 



Quantity 


Magnitude 


Mean Field 


Correlation 


% CI 


Exciton total energy, Ei,± (eo,/io) 

oiexciLoii LOLcii eneigy, -C/2,2 

Total energy of two electrons, Eo,2 (eo) 


1136.3 
335.0 


1138.3 

9977 3 

335.8 


2.0 
i n q 

0.8 


100.1 
i nn 9 

100.1 


Exciton transition energy, Ex, [Eq. (|l])] 
Exciton binding energy, Ax, [Eq. (0)1 

LI i—i 

Biexciton binding energy, Axx, [Eq. (ph] 


1136.3 
46.2 
6.2 


1138.3 
44.1 
-0.6 


2.0 
2.0 
6.8 


100.1 
97.8 
64.5 


1st exciton charging energy, w[ X \ [Eq. (Eh] 
2nd exciton charging energy, W^, [Eq. (Eh] 
1st exciton addition energy, , [Eq. (H)] 


1136.3 
1130.1 

-6.2 


1138.3 
1139.0 
0.6 


2.0 
8.9 
6.8 


100.1 
100.2 
64.5 


1st electron charging energy, fi[ e \ [Eq. (n)] 
2nd electron charging energy, (1% , [Eq. (B)] 
1st electron addition energy, Aj j, [Eq. (pi)] 


147.5 
187.5 
40.0 


147.5 
188.3 
40.8 


0.0 
0.8 
0.8 


100.0 
100.1 
101.4 



rules of Metropolis Monte Carlo integration. The se- 
quence of configurations, Ri, R4+1, . ■ . , samples the den- 
sity |* T (i?)| 2 . 

The real power of QMC is that it can go beyond 
the variational formalism and actually project the true 
ground state energy from an input variational trial func- 
tion, ^t-E^I By weighting the the configuration as it sam- 
ples configuration space, the random walk can identified 
with the imaginary time propagate* exp(— Ht). In this 
diffusion Monte Carlo algorithm,! 5 if 6 ! the random walk in 
configuration space actually samples \&^<I>o where $0 is 
the true ground state wavefunction. The energy expec- 
tation value along the walk E = (ty T \H\$ ) / ($>t\$o) 
is then the true ground state energy of the many-body 
Hamiltonian. That is, even though the true ground state 
wavefunction $0 is never explicitly calculated, its energy 
can be sampled from a random walk. In the remainder of 
the paper, the term QMC will refer to the diffusion Monte 
Carlo algorithm, unless explicitly noted otherwise. 

Applications of QMC to quaafapa e ^ pt s-ha^ie used varia- 
tional QMC, Ell diffusion QMC,BBBNk§EZl and a path- 
integral formulation, related to the diff us i on a lgorithm 
and. based on Feynman path integrals. H'I3e3eII Harju et 
alEil have used both direct diagonalization and VMC to 
calculate the ground state energy of up to 6 electrons in 
a two-dimensional harmonically confined dot. Diffusion 
QMC within the EMA has been used (1) by AustmB 
to calculate the binding energy of excitons inj-a. spherical 
dot as a function of dot radius, (2) by BoltonEa to calcu- 
late the energy of up to 4 electrons in a two-dimensional 
harmonically confined dot ip-the presence of a magnetic 
field, (3) by Shumway et a£J to calculate total energies 
for electron addition to a pyramidal dot, (4) by Pederiva 



et alo to calculate ground and excitation energies for 
up to 13 electrons in a three-dimensional harmonically 
confined dot .and compare to HF and LSDA, and (5) by 
Luczak et a£3 to study energies of up to 20 electrons 
confined to a two-dimensional harmonic potential. Lee 
et aO have used QMC within the EMA to study a pair 
of electrons in a two-dimensional parabolic confining po- 
tential. Path integral QMC has been used by Egger et 
alEj to studied crossover from Fermi liquid to Wigner 
molecule behavior using PIMC within the EMA on up to 
8 electrons in a two-dimensional harmonically confined 
dot, and by Harting et aEJ to calculate the total energy 
of up to 12 electrons in a two-dimensional harmonically 
confined dot. 



2 Configuration Interaction 

In the CI approach, the solutions of the many-body 
Hamiltonian are expanded in terms of Slater determi- 
nants |$) obtained by removing M electrons from the 
valence band and adding TV electrons to the conduction 
band: 

l*> = X] A ( fl l--- h M,ei---e N )\®h 1 ...hM,ei...e N ), 

hi. .Mm ei...e N 

(13) 

where: 

l$fc 1 ...fc M ,e 1 ... ew ) = 4i • • ■ d L c L • ■ -4 N i*o>. (14) 

Here dl ■ ■ ■ dt create holes in the valence-band states 

fti tlM 

hi . . . Hmi while cj • • • cj create electrons in the con- 
duction band states e\ . . . ejv- The Hamiltonian is then 
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diagonalized in the basis of Slater determinants |$) . This 
approach gives access to not only the ground state of the 
system, but also excited states. 

Full CI (FCI) includes all possible determinants from 
a given (finite) set of single particle basis functions, i. e. 
Nh hole orbitals and N e electron orbitals. In the limit of 
an infinite set of basis functions, (Nh, N e ) — > (00,00), 
FCI provides the exact many-body solution, which is 
equivalent to the QMC results. However, most CI ap- 
plications use a small and finite basis set to solve the 
Schrodinger problem. Thus, even including in the CI 
expansion all possible Slater determinants from a finite 
number of single-particle states (FCI) does not provide an 
exact solution, in contrast to QMC. For our calculations 
we also only use a small, finite basis set of bound states, 
denoted (Nh,N e ), therefore ground state total energies 
from FCI will be above the true ground state total en- 
ergy. A useful truncated CI basis is Singles and Doubles 
Configuration Interaction (SDCI), which the set of all 
determinants obtained by exciting at most two particles 
(electrons or holes) from the ground-state (or reference) 
determinant. SDCI is equivalent to FCI for a single ex- 
citon (or two electrons), but is an approximation for two 
or more excitons (or three or more electrons). 

The CI method has been used in the past to solve 
the the many ^hody-. . Schipd inger equation in the EMA, 
approximatiorJjEj'HSEiEl and also tight binding.EHl 
More recently, the CI approach has been used in the con- 
text of the empirical pseudopotential method (EPMlihi. 
single excitonSjEEl electron and hole addition energies, cltB 
and multiexcitonsBj 



III. APPLICATION OF QMC AND CI TO A 
SINGLE-BAND EFFECTIVE-MASS DOT WITH 
FINITE BARRIER 

We first use a simplified single-band EMA model which 
can be treated by both QMC and CI. Our reference 
system is a spherical dot with radius R = 40 A, effec- 
tive masses m e = 0.1 and uih = 0.5, dielectric constant 
e = 12, and barriers AE V = 0.4 eV and AE C = 0.2 eV. 
The energies of the optical and electronic properties of 
this dot are summarized in Table [0]. We have then var- 
ied the radius from to 80 A, while keeping the bar- 
riers fixed. This yields a range of bound electron and 
hole states. The energies of the lowest (i. e. band-edge) 
states eo and ho as a function of dot radius R are shown 
in Fig. ||(a). When the radius R of the dot goes to in- 
finity we have a 3D bulk material called "material I" 
with m e = 0.1, rrih = 0.5, and e = 12. When the ra- 
dius R of the dot goes to zero we have a 3D bulk mate- 
rial called "material II" with m e , m^, and e identical to 
"material I." The band offsets between the two materials 
AE h = 0.2 eV for the valence band and AE e = 0.4 eV 
for the conduction band, so that the band-gap of "mate- 
rial II" is AEh + AE e = 0.6 eV larger than the band-gap 
of "material I." The bulk exciton in both materials is the 
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FIG. 2: Exciton and biexciton binding energy versus dot ra- 
dius as calculated by QMC, for the dot geometry shown in 
the inset. Panel (a) shows the energies of the non-interacting 
electron and hole band edge states. Panel (b) shows the the 
exciton binding energy Ax [Eq. (^)], and biexciton binding 
energy Axx [Eq. (^)] . The bulk exciton Rydberg energy and 
Bohr radius are denoted ao = 7.6 nm and Er = 7.9 meV, 
respectively. Contributions to exciton and biexciton binding 
energy versus dot radius are shown in (c) and (d), respectively. 
Contributions are from: first order perturbation theory (PT), 
self-consistency correction (SC), and correlation (Corr.). 

same, and has a radius ao = 76.2 A, a binding energy 
Eji = 7.873 meV. Both bulk materials have a bound biex- 
citon with the same binding energy, A^J^ = 0.716 meV 
= 0.9Er, (calculated by QMC). In some calculations we 
have varied the barrier energy from AE V = 0.05 eV to 
AE V = 4 eV and AE C = 0.025 eV to AE C = 2 eV, while 
keeping the radius fixed at 40 A. Our model system has 
thus been chosen to roughly capture some properties of 
small SK or colloidal dots, as summarized in Table |n| 
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CI convergence for total energy of an exciton, biexciton, and charged dot 
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FIG. 3: CI convergence of the total energy for three cases: (a) an exciton, (b) a biexciton, and (c) two electrons. All energies are 
measured relative to the center of dot gap. For our CI expansion, we have used single and double substitutions (SDCI) and also 
all possible determinants (FCI). Note that SDCI is equivalent to FCI for cases (a) and (c). SDCI gives a good approximation 
to FCI for case (b), and involves far fewer determinants (see Table III). In all cases our CI expansion captures about half of 
the correlation energy. The correlation energy (and hence CI error) is a very small fraction (< 1%) of the total energy in all 
three cases. 



A Total energies for occupation by an exciton, 
biexciton, and two-electrons 

Figure || shows the total energy for (a) exciton, 
JSi^feJ, hg); (b) biexciton, £2,2(60, h^); and (c) two- 
electrons, £2,o( e o> ^o)- We have decomposed the total 
energies into the three parts listed in Eq. (|§|): first- 
order perturbation theory (_Epx), self-consistent mean 
field (E PT + AE SC ), and the exact QMC result (E tot = 
E PT + AE SC + AE COII ). We then plot the results of CI 
calculations as a function of the number of single-particle 
states (Nh,N e ) used to generate the CI basis set, taking 
either singles and doubles only (SDCI) or all possible de- 
terminants (FCI). The CI energies for one determinant 
are equivalent to the MF result, and the FCI values must 
reach the QMC result in the limit of an infinite basis. The 
total number of CI determinants for M holes and N elec- 
trons occupying Nh hole states and N e electron states is 
C\f h ■ C 2 N N % where C£ = n\/[m\(n - m)!]. The factors 
of 2 are due to the spin-degeneracy of the single particle 
states. Table III lists the actual number of determinants 
for each of the FCI and SDCI data points in Fig. 0. The 
first three lines of Table || give a summary of the role of 
correlation energy and CI convergence in the total energy 
of these three systems. 

In each system, the total energy estimated by first- 
order perturbation theory is above the true ground state 
energy (as required by the variational principle). Self- 



consistency improves upon first-order perturbation the- 
ory, and correlation provides additional improvement. 
For excitons, the self-consistency decreases the energy 
by ~ 1 meV, and correlation gives another ~ 2 meV 
improvement. The total energy, however, is E\ \ — 
1136 meV. So, although our CI only recovers about half 
of the correlation energy, the total energy is only overes- 
timated by about 0.1%. For the case of a biexciton, self- 
consistency also lowers the energy by ~ 2 meV, while cor- 
relation lowers the energy by another ~ 10 meV.— In cal- 
culations on a strain induced dot, Brasken et alEil found 
that SDCI captured &: 90% of the correlation energy for 
multi-excitons, based on comparison to FCI for one to 
four excitons. In our biexciton calculations, we also find 
that SDCI recovers nearly as much correlation energy as 
FCI, but this represents only about about half of the total 
correlation energy. Again, though, correlation represents 
a small part of the total energy of the biexciton, so CI 
(FCI and SDCI) only overestimate the total energy by 
~ 0.2%. For a dot containing two electrons, corrections 
beyond first-order perturbation theory are much smaller, 
~ 1 meV. In fact, for the system calculated here, we 
find only a 0.35 meV decrease in the two-electron system 
with self-consistency, and correlation decreases the total 
energy by about another 0.8 meV. Our CI expansion 
again captures about half this correlation energy, lead- 
ing to negligibly small overestimation of the total energy 
(< 0.1%). 
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TABLE III: Number of determinants used for each of the CI 
calculations shown in Figs. |^ and [|, using only single and dou- 
ble substitutions (SDCI), or all possible slater determinants 
(FCI). Note that SDCI is equivalent to FCI for the case of an 
exciton or two electrons. For FCI the number of CI determi- 
nants for M holes and N electrons occupying Nh hole states 
and N e electron states is C 2 J? h ■ G 2 ?' . 



System 

Exciton (/^e 1 ) 



Biexciton (h 2 e 2 ): 



Two Electrons(/i°e 2 





# of determinants 


(N h ,N e ) 


SDCI 


FCI 


(1,1) 


4 


4 


(4,1) 






(4,4) 


64 


64 


(9,4) 


144 


144 


(10,4) 


160 


160 


(17,4) 


262 


262 


(20,4) 


320 


320 


(1,1) 


1 


1 


(4,1) 


28 


28 


(4,4) 


199 


784 


(9,4) 


564 


4284 


(10,4) 


649 


5320 


(17,4) 


1356 


15708 


(20,4) 


1719 


21840 


(0,1) 


4 


4 


(0,4) 


28 


28 


(0,5) 


153 


153 


(0,10) 


190 


190 



B Exciton and biexciton transition and binding 
energies 

Measured quantities such as the exciton and biexci- 
ton binding energies represent differences between total 
energies. Even if the mean- field contributions dominate 
total energies, the mean-field contributions to differences 
of total energies may have significant contributions from 
correlation. Lines 4-6 of Table [n] summarize the role 
of correlation and CI convergence for the exciton tran- 
sition energy, Ex [Eq. Q]; exciton binding energy, Ax 
[Eq. (g)]; and the biexciton binding energy, Axx [Eq. 
(|J)]. Correlation is only a small part (2 meV) of the 
exciton transition energy Ex = 1136.3 meV. So, even 
though our underconverged CI fails to capture all the cor- 
relation energy, Ex is only overestimated by 0.1%. The 
same 2 meV of correlation energy is a much larger com- 
ponent of the exciton binding energy, = 46.2 meV, 
so errors due to underconvergence of CI are more sig- 
nificant, and CI underestimates Ax by more than 2%. 
The biexciton binding energy Axx = 6.2 meV is due 
entirely to 6.8 meV of correlation energy, so CI under- 
convergence is much more serious. Our CI calculation of 
biexciton binding is only 65% of the exact QMC result. 

In Fig. ^ we show the results of first-order perturbation 



theory (£pt), self-consistent mean field (Ept + AEsc), 
the exact QMC result (E tot = E PT + AE SC + AE co „), 
and CI convergence vs. basis size for (a) the exciton tran- 
sition energy and (b) the biexciton binding energy. For 
the exciton transition energy, Fig. 0(a), increasing the CI 
basis does improve the calculated energy, but it is only 
a difference of ~ 2 meV out of a much larger exciton 
transition energy of 1.136 eV. On the other hand, the 
CI correction is essential to even approximate the biex- 
citon binding energy, shown in Fig. [|(b). Note that the 
improvement of the biexciton binding with CI basis size 
is not monotonic. This is because the biexciton binding 
is a difference of one- and two-exciton energies. As the 
basis is increased, the relative improvement in the one- 
and two-exciton total energies varies, thus the calculated 
biexciton binding energy can actually decrease when the 
CI basis is improved. We also show the results of SDCI 
in Fig. |(b). 



1 Dependence on dot size 

We have varied the dot radius from R = to R = 80 A, 
all in the strongly confined regime, R . ao — 76.2 A. 
Figure |2|(b) shows the exciton and biexciton binding en- 
ergies as calculated by QMC. Figures |(c) and |(d) de- 
compose the contributions to the exciton and biexciton 
binding into (1) first order perturbation theory, (2) self- 
consistency corrections, and (3) correlation corrections, 
as in Eq. (||). 

The small R limit is the energy of a bulk-II material, 
and all excitonic binding energy is from correlation. As 
the radius of the dot increases, the bulk-II exciton binds 
to the dot, the exciton binding energy is enhanced, and 
most of the binding energy comes from perturbation the- 
ory. The maximum in the binding energy occurs when 
the electron and hole are both individually bound to the 
dot, but the radius is small, so that the direct Coulomb 
interaction (from first order perturbation theory) is the 
strongest. The exciton binding energy exhibits a clear 
peak at around R ~ 40, A, in similarity with previous 
calculations by Austin.E3 As the dot becomes larger, the 
direct Coulomb interaction from perturbation theory de- 
creases, causing a decrease in the exciton binding en- 
ergy. Finally, as the dot becomes comparable in size to 
the bulk-I exciton radius, correlation begins to have sig- 
nificant contributions to exciton binding. In the limit 
R ^ ao (not shown) , the binding energy becomes that of 
a bulk-I exciton. 

The biexciton binding energy is greatly enhanced in 
a quantum dot, except for the case of a very small dot 
with only a single weakly bound exciton. We find that 
the biexciton binding energy is remarkably insensitive to 
dot radius, having a value Axx between 5.1 meV and 
6.2 meV (0.7 E R to 0.9 E R ) for dots with radii R between 
2 nm and 8 nm (0.3 ao and 1.1 ao). This is in contrast 
the exciton binding energy, A^, which exhibits a clear 



9 



CI convergence for exciton transition, biexciton binding, and electron charging energies 
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FIG. 4: CI convergence of: (a) exciton transition energy, Ex [Eq. (Eh]; (b) biexciton binding energy, Axx [Eq. (^)]; and (c) 
second electron addition energy, /j,^ [Eq. (0)]. For our CI expansion, we have used single and double substitutions (SDCI) 
and also all possible determinants (FCI). Note that SDCI is equivalent to FCI for cases (a) and (c). SDCI gives a good 
approximation to FCI for case (b), and involves far fewer determinants (see Table |[Il| ). In all cases our CI expansion captures 
about half of the correlation energy. The correlation energy (and hence CI error) is a significant fraction of the the total energy 
only for case (b), biexciton binding. 



peak at small dot radius. The size range 10 . R . 18 A 
has a negative biexciton binding. Physically, these are 
small dots that can weakly bind two excitons, but with a 
higher total energy than separating the two excitons on 
two non-interacting, identical dots. We see from Fig. ||(d) 
that the biexciton binding energy is almost entirely due 
to correlation, as noted before. 

2 Dependence on barrier height 

To study the effect of finite confining barriers on ex- 
citon and biexciton binding energies, we have varied the 
dot barriers from zero to infinity. In all calculations we 
have kept AE e /AEh = 2 and used a radius of 40 A. In 
Fig. ||(b) we plot the binding energies of excitons and 
biexcitons calculated with QMC as a function of barrier 
height. The 40 A dot is able to bind an electron once 
AE e &: 30 meV, and binds a hole once AEh &; 5 meV. 
Unlike the behavior seen with varying the dot radius, 
increasing the confining potential leads to a monotonic 
increase in exciton and biexciton binding energies. For 
zero barrier potential, the exciton has the bulk-I exciton 
binding energy, Ax = EjP = 7.9 meV. As the barrier 
potential is increased enough to bind both electrons and 
holes, the exciton binding increases rapidly. The binding 
energy reaches a maximum of Ax = 55 meV = 7Er 
for infinite barriers. Similarly, the biexciton binding 



energy starts from the bulk biexciton binding energy 
Axx = 0.7 meV = O.IEr and increases to a maximum 
of Axx = 7.2 meV = 1.0-Er for infinite barriers. Fig- 
ures |^(c) and |](d) show the contributions of perturbation 
theory, self-consistency correction, and correlation to the 
exciton and biexciton binding energy. Except for very 
weakly confined dots, the exciton is very well described 
by first-order perturbation theory. For weak confine- 
ment, the electron is unable to bind, but self-consistent 
interaction with the hole is able to bind the electron, so 
that the exciton binding energy is almost entirely due 
to self-consistency. For the weakest confinement, nei- 
ther the electron nor the hole is bound, and the excitonic 
binding is entirely due to correlation. Again, biexciton 
binding is due entirely to correlation. 

C Multi-exciton energies 

Figure || shows mean- field and exact (QMC) results for 
the multi-exciton charging energies Wjv [Eq. (Q)], and 

the multi-exciton additions energies, A^-]^ +1 [Eq. (^)]. 
The most prominent feature is the jump in the charg- 
ing energy for W3, which also appears as a peak in the 

(x) 

addition energy A;, 3 . This "shell effect" arises because 
only the first two excitons can occupy the lowest energy 
eo and ho states. Starting with the third exciton, Pauli 
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FIG. 5: Exciton and biexciton binding energy (including cor- 
relation) as calculated by QMC as a function of versus bar- 
rier energy, with the constraint AE e / AEh = 2, for the dot 
geometry shown in the inset. Panel (a) shows the single- 
particle energies of the non-interacting electron and hole band 
edge states. Panel (b) shows the the exciton binding energy 
Ax [Eq. (^)] and biexciton binding energy Axx [Eq. (0)]. 
The bulk exciton Rydberg energy is denoted Er — 7.9 meV. 
Contributions to exciton and biexciton binding energy versus 
barrier energy are shown (c) and (d), respectively. Contri- 
butions are from: first order perturbation theory (PT), self- 
consistency correction (SC), and correlation (Corr.). 



exclusion requires the addition excitons to start filling 
the next energy shell, e\h\ through 63/13. This is a fea- 
ture of the single particle model, and does not require 
any treatment of correlation. Correlation is necessary to 
describe the decrease in charging energy for the second 
exciton, W2 < W\, or equivalently the negative value of 

(-X) 

the first exciton addition energy A\ 2 ' = —6.2 meV. This 
is the positive biexciton binding energy Axx = 6.2 meV, 
discussed earlier. As shown in lines 7-9 of Table |1 the 
correlation contribution for the second charging energy 
W2 is 8.9 meV, considerably larger than the 2.0 meV 
for Wi . Our CI only captures about half the correlation 
energy, so it slightly overestimates the exciton charing 
energies, and considerably underestimates the negative 
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FIG. 6: (a) Exciton charging [Eq. (||)], and (b) addition ener- 
gies [Eq. (rj)], for the dot shown in the inset. Because excitons 
are neutral, it is energetically favorable for a dot to hold many 
excitons. 



value of A^'. 



D Electron loading energies 

Figure shows mean- field and exact (QMC) results for 
electron charging energies /Ltjv, [Eq. (||)], and the electron 

additions energies, A^- ) Ar+1 [Eq. (0)]. Because electrons 
are charged, Coulomb repulsion quickly limits the num- 
ber of electrons that can be loaded into the dot. For 
our model, shown in the inset to Fig. [7], it is only en- 
ergetically favorable to add four electrons; beyond this, 
electrons would rather escape into the barrier material 
conduction band, shown as a dashed horizontal line in 
Fig. 0(a). The re is a peak in the electron addition en- 
ergy A\ 3 in Fig. 0(b). This is due the filling of the Bq 
state by a spin-up and spin down electron (another "shell 
effect"). Both QMC and MF capture this single particle 
effect. As shown in Fig [|(c), our CI expansion recov- 
ers about half of the correlation energy for two electrons. 
However, the correlation energy in a two-electron dot is 
only about 1 meV, so CI errors are a negligible 0.5 meV. 
The small value of correlation and the good agreement 
of our CI calculations for dot charging are summarized 
in last three lines of Table |l] 
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FIG. 7: (a) Electron charging energies [Eq. (^)], and (b) ad- 
dition energies [Eq. (Q)], for the dot shown in the inset. This 
dot can only hold up to four electrons, due to Coulomb re- 
pulsion. The conduction band minimum energy of the barrier 
material, AE e = 400 meV is shown in panel (a). 



IV. APPLICATION OF CI TO A MULTI-BAND 
DOT DESCRIBED VIA PLANE- WAVE 
PSEUDOPOTENTIALS 

QMC calculations are currently limited to either 
small systems containing up to a few hundreds of 



electrons or to highly simplified model Hamil- 
tonians (such as the EM A). A more accurate descrip- 
tion of the electronic structure (Fig. [l]) of semiconductor 
quantum dots can be obtained using the pseudopotential 
approach.H Unfortunately, QMC methods are presently 
unable to deal with the large number of electrons of a typ- 
ical quantum dot, and CI is the only viable approach to 
treat correlation effects in large quantum dots described 
by atomistic pseudopotentials. In addition, the diagonal- 
ization of the CI Hamiltonian gives access to the excited 
states (unavailable in ground-state QMC calculations) as 
well as the ground state of the electronic system, thus en- 
abling the calculation of the optical spectrum of quantum 
dots. 

In order to illustrate the capabilities of the CI ap- 
proach combined with a pseudopotential description of 
the electronic structure, we consider a nearly spherical 
CdSe quantum dot having the wurtzite lattice structure 
and a diameter of 38.5 A. The surface dangling bonds are 
fully passivated using ligand-like atoms. eZI This quantum 
dot is representative of CdSe nanocrystals grown by col- 
loidal chemistry methods. 

We consider here only low-energy excitations of the 
electronic system, which are obtained by promoting elec- 
trons from states near the top of the valence band to 
states near the bottom of the conduction band. The 
band-edge solutions of Eq. ( |l0| ) can be efficiently ob- 
tained using the folded spectrum method,Hoo which 
allows one to calculate selected eigenstates of the 
Schrodinger equation with a computational cost that 
scales only linearly with the size of the system. In this 
approach, Eq. ( |l0| ) is replaced by the folded-spectrum 
equation 



r - i 2 

-V 2 +V ps (r)+V NL -£ref 

where e re / is an arbitrary reference energy. The lowest 
energy eigenstate of Eq. (|l^) coincides with the solu- 
tion of the Schrodinger equation [Eq. (]To|) ] whose en- 
ergy is closest to the reference energy e re /. Therefore, by 
choosing the reference energy in the band gap, the band 
edge states can be obtained by minimizing the functional 

The solution of Eq. ( pTq ) is performed by expanding 
the wave functions tf>i(r, a) in a plane- wave basis set. 
To this purpose, the total pseudopotential V ps (r) is de- 
fined in a periodically repeated supercell f2 containing 
the quantum dot and a portion of the surrounding ma- 
terial. The supercell £1 is sufficiently large to ensure that 
the solutions of Eq. (Jlf^) are converged within 1 meV. 
The single-particle wave functions can then be expanded 
as ipi(r, a) — J2g c «(G,cr) exp(iG • r), where the sum 
runs over the reciprocal lattice vectors G of the supercell 
Q. The energy cutoff of the plane-wave expansion is the 



>i(r,(T) = (e°i -e ref f Vi(r,<r), (15) 

same used to fit the bulk electronic structure, to ensure 
that the band-structure consistently approaches the bulk 
limit. The minimization of the functional A[ip] is carried 
out in the plane-wave basis set using a preconditioned 
conjugate-gradients algorithm. 

In the next step we construct a set of Slater determi- 
nants \$h 1 ---h N ,e 1 ---e N ) [see Eq. (O)] obtained by creat- 
ing N holes in the valence band and N electrons in the 
conduction band, and diagonalize the CI Hamiltonian in 
this basis set. Using the CI approach, we have calculated 
the multiexciton spectrum of a CdSe dot. We consider 
here up to three excitons and we use a CI basis set of 
480 configurations for the single exciton, 43890 config- 
urations for the biexciton, and 20384 configurations for 
the tri-exciton. All the relevant interactions (including 
electron-hole exchange) are included in the CI calcula- 
tions. We assume that when an A^-exciton is created in 
the quantum dot, it relaxes non-radiatively to the ground 
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FIG. 8: Schematic illustration of the leading contributions 
to peaks (A3, B3, A' 2 , A2, Ai) appearing in Fig. Solid hor- 
izontal lines are energies of N = to N = 4 excitons, with 
dashed lines indicating states that do not participate in dipole 
transitions. 

state before decaying radiatively into an (JV— l)-exciton. 

The calculated multiplet levels are shown in Fig. |^ and 
the emission spectrum is shown in Fig. ^. The three 
panels of Fig. ^ correspond to the recombination of (a) 
a tri-exciton into a biexciton (3 — > 2), (b) a biexciton 
into a single-exciton (2 — > 1), and (c) a single-exciton 
into the ground state (1 — > 0), respectively. We assume 
that the low-energy states of the 7V-exciton are thermally 
populated (kT — 5 meV) before recombination. We see 
from Fig. ^ that: 

(i) The single-exciton recombination spectrum, 
Fig. ||(a), shows a single peak (Ai) centered at 2.154 eV. 
It is well knownO that in CdSe nanocrystals the 
electron-hole exchange interaction splits the lowest- 
energy excitonic state (/iQ,eJ) into two doublets, having 
total angular momentum F — 2 and F = 1 respectively 
(see Fig. || ). The lower-energy doublet (F — 2) is 
optically forbidden, while the higher-energy doublet 
(F = 1) is optically allowed. We find an energy sep- 



aration of ~ 5 meV between the two doublets. The 
emission peak A\ observed in Fig. ^ comes from the 
recombination of the higher-energy doublet, which is 
thermally populated. This explains the relatively weak 
intensity of the single-exciton peak. 

(ii) The biexciton recombination spectrum, Fig. ^(b), 
shows a strong peak (A2) centered at 2.140 eV. This 
peak originates from the recombination of a biexciton 
in the ground state (/ig,e§) into a single exciton in the 
F = 1 state. The weak shoulder to the red of the main 
peak (A'2) is due to the recombination of a thermally oc- 
cupied higher-energy biexciton state in the configuration 
e o)- Note that several transitions from the biex- 
citon ground state to single-exciton excited states are in 
principle possible, but have very weak oscillator strength. 
These transitions would occur to the red of the fundamen- 
tal transition. The calculated biexciton binding energy 
is 2Ex — Exx ~ 4 meV. This value is probably un- 
derestimated due to the under-convergence of the CI ex- 
pansion. Interestingly, the "apparent" biexciton binding 
energy, i.e. the red-shift of the main biexciton peak A 2 
with respect to the single-exciton peak A\, is ~ 14 meV 
(not 4 meV). The reason is that the biexciton recombi- 



15 
10 
5 




(a) 



A 3 



3^2 



: 20 




2.2 2.3 2.4 

Energy (eV) 



2.6 



FIG. 9: Exciton transition energies for a CdSe dot for (a) 
decay from three to two excitons, (b) decay from two to one 
excitons, and (c) decay of a single exciton. The intensity scale 
is different in each of the three panels, and weak transitions 
between peaks A3 and B3 in (a) have been magnified by x 20. 
Grey vertical lines indicate all calculated transition energies, 
and solid black lines in the Gaussian broadened transitions 
weighted by calculated dipole transition strengths. 
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nation takes the quantum dot in the F = 1 excited state, 
rather than the F — 2 ground state (see Fig. ||). Thus we 
haveE(Ai)-E(A 2 ) = (E^ 1 - E . )-{E XX ~ E^ 1 ) = 
A X x + 2 (£f =1 - E^= 2 ) =4 + 2x5 meV = 14 meV. 

(iii) In the case of three excitons we find that the 
ground state wave function originates primarily from the 
non-Aufbau configuration }v^h\\e\e\. In fact, the third 
hole prefers to occupy the p-like h% state rather than the 
s-like hi state, due to reduced Coulomb repulsion with 
the remaining two holes. Two main transitions are pos- 
sible from the three-exciton ground state: the eo — > h a 
recombination, which leaves the system into the excited 
biexciton configuration /ij/i3,;eje}, leads to peak A3 lo- 
cated at 2.188 eV. The ei — > hi recombination, which 
takes the system into the ground-state biexciton con- 
figuration fto! e O! ^ s responsible for peak B3 centered at 
2.497 eV. Note that the B3 transition originates from 
an exchange-split tri-exciton state (see Fig. ||) which is 
thermally populated, hence the relatively weak oscillator 
strength of the B 3 transition. 

Note that a calculation considering only ground-state 
to ground-state transitions would miss most of the peaks 
observed in Fig. ^. The capability of the CI expansion 
to access excited states, coupled with the possibility of 
using a multi-band pseudopotential Hamiltonian for the 
calculation of the single-particle energies and wave func- 
tions, makes it the method of choice for calculating ex- 
cited states of semiconductor quantum dots. 

V. CONCLUSION 

We have studied the effects of correlation on a simpli- 
fied, single-band model dot using both QMC and CI, and 
have studied correlation in the multi-exciton PL spectra 
of a realistically modeled CdSe dot using CI. Our results 
for the simplified, single band model are summarized in 
Table ||. We find the following results for our model: (1) 
total energies for an exciton, biexciton, and two electrons 
are dominated by mean field effects, so that correlation 
energies and CI convergence errors are less than 1% [see 
Fig. H]; (2) typical exciton transition energies, which are 
~ 1 eV, can be calculated to closer than 1% by perturba- 
tion theory, with only a ~ 2 meV correlation correction 
[see Fig. Q(a)]; (3) typical exciton binding energies are 
~ 46 meV, with only 2 meV from correlation, and our 
CI captures roughly half of the correlation to give ex- 
citon binding energies that are nearly 98% of the exact 
QMC value; (4) typical biexciton binding energies are 
positive ~ 6 meV, almost entirely due to correlation en- 
ergy, and our CI only recovers about 65% of the exact 
QMC value [see Fig. [|(b)]; (5) exciton charging energies 
are ~ 1130 meV and well described by CI, while exci- 
ton addition energies can be due entirely to correlation, 
in which case our CI is only qualitatively correct; and 
(6) typical electron charging energies are ~ 150 meV, of 
which correlation contributes very little (~ 1 meV) , like- 



wise, electron addition energies are ~ 40 meV with very 
little correlation contribution, so that CI is accurate to 
about 1-2% for electron addition energies. 

Although QMC is a good method for testing conver- 
gence of CI on a simplified, single band model, only CI 
may be used on our more realistic model of CdSe. Our 
multi-band pseudopotential model capture the correct 
symmetries and electronic structure of the dots, lead- 
ing to qualitatively different predictions than single-band 
models. For example, the multiplet structure presented 
in Fig. ^ requires a multi-band description of the sin- 
gle particle levels. Some of the details of our realistic 
CdSe calculation that are missing from our single-band 
CI model are: (1) different degeneracies of the single- 
particle hole levels due to a multi-band description of the 
valence band states, (2) electron-hole exchange splitting 
of 5 meV in the ground state (/ij,ej) exciton, (4) the 
existence of many weak transitions that are symmetry 
forbidden in single band models, An additional benefit 
of CI is that it gives excited state energies necessary to 
identify some of the peaks that appear in single-dot pho- 
toluminescence spectra. 

We conclude that correlation effects are important to 
some quantities, such as exciton binding and exciton ad- 
dition energies, and essential to calculate positive bind- 
ing energies. QMC methods are well-suited for sim- 
ple, single-band models. Applications to realistic mod- 
els which capture the proper symmetries and electronic 
structure of quantum dots are currently restricted to 
CI methods. We find that CI calculations including all 
bound states are accurate to better than 3% for many 
measurable properties, as listed in Table |l[ Even for 
biexciton binding, which is dominated by correlation, 
our CI calculations are qualitatively correct, capturing 
about 65% of the QMC prediction for a simplified mod- 
els. Therefore we conclude that realistic multi-band mod- 
els combined with perturbation theory and a judicious 
use of CI for correlation corrections is a computational 
approach well-suited to realistic modeling of interacting 
electrons and holes in SK and colloidal semiconductor 
quantum dots. 
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